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The neutrino experiment KARMEN at the beam stop neutrino source ISIS investigates the oscillation channel 
Vfi^ in the appearance mode by looking for p {ve ) n reactions. An analysis of data collected from February 
1997 through March 2000 with the KARMEN2 experimental setup reveals 11 candidate events in good agreement 
with the background expectation of 12.3 ± 0.6 events. Hence, there is no indication of an oscillation signal. A 
maximum likelihood analysis of the data leads to an upper limit (at 90% confidence level) for the mixing angle 
of sin^(2e) < 1.3 • 10"^ at large Am^ and Am^ < 0.049 eVVc" for sin^(2e) = 1. 

The anomaly in the time spectrum of events induced by and i/p seen in the KARMEN 1 data could not be 
confirmed with the KARMEN2 data. 



1. THE EXPERIMENT 

The KARMEN experiment is performed at the 
neutron spallation facility ISIS of the Rutherford 
Appleton Laboratory. Neutrinos are produced by 
stopping 800 MeV protons in a massive beam stop 
target, thereby producing pions. The 7r~ are ab- 
sorbed by the target nuclei whereas the n'^ de- 
cay at rest (DAR). Muon neutrinos therefore 
emerge from the decay 7r+— > i^^,. The low 

momentum /i+ are also stopped within the mas- 
sive target and decay via /^^ — > e"*" -|- i/g + • 
Because of this 7r+-/x+-decay chain at rest ISIS 
represents a z^-source with identical intensities 
for f^, Ue and emitted isotropically ($y = 
6.37 ■ lO^^v/s per flavor for a proton beam cur- 
rent Ip — 200 /xA). There is a minor fraction of 
TT~ decaying in flight (DIE) with the following 
fj,~ DAR in the target station which again is sup- 
pressed by muon capture of the high Z material 
of the spallation target. This decay chain leads to 
a very small contamination of v^jv^. < 6.2 • 10^^ 
|2|, which is further reduced by software cuts. 

The energy spectra of the i/'s are well defined 
due to the DAR of both the 7r+ and /i"^ (Fig- 
ure |l|a). The v^^^ from tt"*" -decay are monoen- 
ergetic with E(i/^)=29.8 MeV, the continuous en- 
ergy distributions of Vf. and i^^^ up to 52.8 MeV 
can be calculated using the V-A theory. Two 
parabolic proton pulses of 100 ns base width and 



a gap of 225 ns are produced with a repetition fre- 
quency of 50 Hz. The different lifetimes of pions 
(r = 26ns) and muons (t — 2.2 p£) allow a clear 
separation in time of the t/^-burst (Figure |l|b) 
from the following j/g's and i^^'s (Figure 0c). The 
accelerator's duty cycle of 10~^ allows effective 
suppression of cosmic induced background. 

The neutrinos are detected in a rectangular 
tank filled with 56 1 of a liquid scintillator ||] . The 
central scintillation calorimeter is segmented into 
512 optically individual modules. The event posi- 
tion is determined by the hit module and the time 
difference of the PM signals at each end of this 
module. Gd2 03 coated paper within the mod- 
ule walls provides efficient detection of thermal 
neutrons due to the very high capture cross sec- 
tion of the Gd(n,7) reaction (a k. 49000 barn) 
in addition to the p ( n,7 ) d capture. The KAR- 
MEN electronics is synchronized to the ISIS pro- 
ton pulses to an accuracy of better than ±2ns, 
so that the time structure of the neutrinos can be 
exploited in full detail. 

A massive blockhouse of 7000 1 of steel in com- 
bination with a system of two layers of active veto 
counters provides shielding against beam corre- 
lated spallation neutron background, suppression 
of the hadronic component of cosmic radiation 
as well as reduction of the flux of cosmic muons. 
In 1996 an additional third veto counter system 
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Figure 1. Neutrino energy spectra (a) and pro- 
duction times of z^^ (b) and I'eji^n (c) at ISIS. 



with a total area of 300 was installed within 
the 3 m thick roof and the 2-3 m thick walls of 
the iron shielding ji) (KARMEN2 experimental 
configuration). By detecting muons which pass 
through the steel at a distance of less than a me- 
ter from the main detector and therefore vetoing 
the successive energetic neutrons from muon deep 
inelastic scattering, the main background for the 
P^— >i/e oscillation search could be reduced by a 
factor of 40 compared to the KARMENl setup. 



2. THE Pf,^ve OSCILLATION SEARCH 

The signature for the detection of is a spa- 
tially correlated delayed coincidence of positrons 
from p ( f'e , e^ ) n with energies up to = 
Ep^-Q = 52.8 - 1.8 = 51.0 MeV and 7 emis- 
sion of either of the two neutron capture pro- 
cesses p(n,7)d with one 7 of £^(7) = 2.2 MeV 
or Gd ( n,7 ) with 3 7~-quanta on average and a 
sum energy of ^ E{'y) — 8 MeV. The positrons 
are expected in a time window of several /iS af- 
ter beam-on-target with a 2.2 fj,s exponential de- 
crease due to the /i+ decay. The time difference 
between the e"*" and the capture 7 is given by the 
thermalization, diffusion and capture of neutrons, 
T„ « 110 ^s. 

The raw data investigated for this oscillation 
search were recorded in the measuring period of 
February 1997 through March 2000 which corres- 
ponds to 7160 C protons on target. A positron 
candidate is accepted only if there is no previous 
activity in the central detector nor in the two in- 
nermost veto counters up to 24 /xs. The required 
cuts in energy and time for the prompt (p) event 
are: 0.6 < tp < 10.6 ^s, 16.0 < Ep < 50.0 MeV. 
The cuts on the delayed expected neutron event 
are as follows: 5.0 < td-tp < 300 ^s, Ed < 8MeV 
and a volume of 1.3 m^ for the spatial coincidence. 
Applying all cuts, the total background expecta- 
tion amounts to 12.3 ± 0.6 sequences where the 
individual background sources are events induced 
by cosmic muons, ^^C {i^e ,e~ ) ^"^^g.s. sequences, 
j^-induced accidental coincidences and (e"'",n) se- 
quences from the intrinsic ISIS De contamination. 

In Table |l] the individual contributions of the 
above described background sources are summa- 
rized. The last row shows the expectation of 
(e+ ,n) sequences from oscillations assuming max- 
imal mixing and Am^ > lOOeV^/c*. The back- 
ground components are determined precisely dur- 
ing the normal measurements: The v induced 
backgrounds are measured with KARMEN in dif- 
ferent energy and delayed time windows, the cos- 
mic background is measured with high statistics 
in the long pre-beam time window. Only the 
small De contribution of the intrinsic source con- 
tamination has to be simulated. In addition, the 
capability of the KARMEN experiment to iden- 
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Figure 2. Spectra of the 11 candidate se- 
quences after applying all cuts. Also shown 
are the background contributions (from bot- 
tom to top): cosmic background, (e~,e"'") from 



12 



C ( I'e , e ) Ng.s., V induced random coinci- 
dences and ISIS Vp contamination. 



Table 1 

Expected sequences from different background 
components within the evaluation cuts specified 
in the text. Last row: oscillation expectation for 
maximal mixing. 



background contribution 


events 


(e-,e+) from ^^C ( i^e , e" ) ^"Ng.,. 


3.9±0.5 


V induced random coincidences 


3.5±0.3 


ISIS Ve contamination 


1.7±0.2 


cosmic induced sequences 


3.2±0.2 


total background 


12.3±0.6 


Ve signal for sin" (20) = 1 


2442±269 



tify neutrino induced events is constantly moni- 
tored by measuring neutrino-nucleus interactions 
on ^^C via neutral and charged current reactions. 

Analysing the data results in 11 sequential 
events which satisfy all conditions (see Figure 
|). This number is in good agreement with 
the total background expectation. Applying a 
Bayesian renormalisation of the physically al- 
lowed region to the experimental result near a 
boundary {N{osc) > 0, 11 measured events with 
12.3 ± 0.6 background events expected), an up- 
per limit of N{osc) < 6.3 at 90% CL can be ex- 
tracted. An identical procedure had been applied 
to the initial KARMEN2 data, when the evalua- 
tion of a potential oscillation signal was based on 
a pure counting experiment due to the very small 
statistics of the background j|] . 

However, with more data and spectral infor- 
mation, a much better evaluation method is a 
maximum likelihood analysis. Such a likelihood 
analysis to extract a possible v^^Vp signal from 
these 11 sequences makes use of the precise spec- 
tral knowledge of all background sources and a 
detailed MC description of the oscillation signa- 
ture in the detector. The likelihood function L 
defined as 

N 
i=l 

is optimised with respect to the free parame- 
ters Am^ and sin'^(20). The probability den- 
sity function / is calculated for each of the 
N event sequences from the parameters x = 
{-^prompt ; ^de layed, tprompt, At, Af) where At and 
Ar denote the delayed spatial coincidence. The 
analysis results in a best fit value of oscillation 
events N{osc) = within the physically allowed 
range of parameters. As shown in Figure ^ an 
upper limit of 3.8 and 3.1 oscillation events for 
Am^ < leVVc** and Am^ > 20eVVc'*, respec- 
tively, can be extracted at 90% CL based on a 
complete frequentist approach as suggested by . 

The hmits at 90% CL and 95% CL can be com- 
pared with the number of oscillation events ex- 
pected from the latest results of the LSND ex- 
periment 0. In Figure ||, the calculated LSND 
signal strength within the KARMEN detector is 
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Figure 3. KARMEN exclusion limits at 90% CL 
and 95% CL in comparison with the range of ex- 
pected oscillation events deduced from the LSND 
evidence 0. The darker (lighter) region corres- 
ponds to logarithmic likelihood values less than 
2.3 (4.6) units below the global maximum. 



shown. For large Am^, even at 95% CL, all of 
the favored LSND signal range can be excluded 
by KARMEN2. At lower values of Am^, a small 
fraction of the signal strength of LSND cannot be 
ruled out by KARMEN2. Note that at high Am^ 
the band of expected events and the exclusion 
curve are almost parallel. With decreasing val- 
ues of Am^ the expected signal strength becomes 
smaller. This may be attributed, in part, to the 
larger source distance for the LSND experiment 
{{(Ilsnd) « 30 m compared to {(Ikarmen) « 
17.7 m) but may also reflect a smaller signal 
strength: If the LSND data evaluation yields an 
oscillation probability P(p^^ 9^) which decreases 
with smaller values of Am^ -instead of remaining 
constant- this would have the same effect. 

Assuming maximal mixing (sin2(2e) = 1), 
2442 ±269 (e+,n) sequences from oscillations with 
large Am^ were expected. Together with the 



limit of 3.1 events deduced in the unified frequen- 
tist approach, this leads to an upper limit on the 
mixing amplitude of 



sin2(2e) < 1.3-10 



-3 



(90% CL) 



for Ato^ > lOOeV^/c''. The corrsponding exclu- 
sion curve in (Am^,sin^(20)) is given in Fig. 0. 
Also shown are limits from other experiments M- 
pT| as well as the favored regions from LSND 
based on a complete re-analysis of the entire 1993- 
98 data set 0. Again, at high Am'^, KARMEN 
excludes the region favored by LSND. At low 
Am^, KARMEN leaves some statistical space, 
but the reactor experiments at Bugey and Chooz 
add stringent limits from the 9^ disappearance 
search. Any exclusion curve or favored region is 
a compactiflcation of a complex statistical infor- 
mation, it reflects only a contour of the twodi- 
mensional likelihood functions which have non- 
trivial properties like multiple side maxima and 
non-parabolic shapes. Therefore, one needs a 
quantitative statistical analysis of both LSND 
and KARMEN based on the detailed event-by- 
event information as demonstrated in [ p^ for pre- 
liminary data sets to deduce correct statements 
of compatibility or disagreement in terms of fre- 
quentist confidence regions. 

Taking the LSND region published first ^ 



which corresponds to a stronger 



signal. 



there is a much stronger disagreement in the ex- 
periments' outcome: The actual 90% CL KAR- 
MEN2 limit excludes the complete parameter 
range favored by the 1993-95 LSND data. 

One way of estimating the sensitivity of an ex- 
periment is to determine the average limit on the 
oscillation parameter simulating a large number 
of experimental outcomes with the actual level 
of background events, but no oscillation signal. 
These samples are subsequently analysed with 
the same maximum likelihood analysis used for 
the real data set. The actual KARMEN limit is 
slightly better than its sensitivity of sin^(20) = 
1.8 • 10~^ for large Arn^ and almost identical at 
lower Am^ (see dashed line in Fig. 4). Compared 
to the earlier results of KARMEN |5|, this corres- 
ponds to an improvement of the sensitivity by a 
factor of « 2.5. 




Figure 4. KARMEN2 exclusion limit and sensitivity (dashed line) at 90% CL compared to other exper- 
iments: BNL E776 H, CCFR §, BUGEY |§, CHOOZ ^ and the evidence for Df,^i?e oscillations 
reported by LSND 



KARMEN will continue to take data for the 
I'tJ.^ search until the end of April 2001 with 
an anticipated sensitivity of sin^(20) = 1.3 • 10^"^ 
for large Am^. Meanwhile, the upcoming BooNE 
experiment at Fermilab is under construction. Its 
sensitivity is expected to improve the final KAR- 
MEN sensitivity by another factor of 2 in the 
appearance mode I'p,—^ Ve with different systema- 
tics than KARMEN and LSND. 

3. THE TIME ANOMALY 

In the time window of 0.6-10.6 /iS after beam- 
on-target, reactions on -"^^C, "'^'^C and ^^Fe induced 
by Ve and are expected to reflect the 2.2 jjLS life 
time of the muons decaying in the target super- 
imposed on a flat, cosmic induced background. 
However, in the KARMENl data there was a 
bump-like distortion of the spectrum at about 



3.1-4.1 /xs after beam-on-target first reported in 
IH]. Taking the total KARMENl data set of so- 
called 'single prong' events (unaccompanied en- 
ergy deposits without delayed coincidences) col- 
lected from July 1990 through December 1995, 
Nx = 89 ± 24 events were observed in excess 
to the neutrino induced events and the cosmic 
background. Detailed analyses and special mea- 
surements did not lead to the identification of po- 
tential background sources(e.g. beam correlated 
neutrons from the spallation source, electronic ef- 
fects, detector inefficiencies or after-pulses of the 
ISIS proton synchrotron) as origin of this excess. 
Therefore, a working hypothesis of a particle pro- 
duction was considered to get a consistent de- 
scription of the origin of these events. 

The ansatz consisted of a rare non-SM pion 
decay mode 7r+^ AT where X denotes a 



massive (mx ~ 33.9 MeV), neutral, weakly inter- 
acting, unstable particle with decay products de- 
tectable via electromagnetic interactions. A pos- 
sible branching ratio as small as 10~^^ relative to 
the standard decay 7r+^ could have ex- 

plained the observed excess without contradiction 
from other experiments. In this framework, the 
X particle would be produced in the ISIS target 
within the pion decay double pulse structure (see 
Fig. I^a), then slowly moving {vx ~ 5m/^s) to- 
wards the KARMEN detector at a mean distance 
of (d) = 17.7 m. This would result in a specific 
double pulse structure in time and position within 
the 3.5 m long scintillation tank. The consistency 
of this hypothesis was tested by a maximum like- 
lihood analysis of individual event times and lo- 
cations . The free parameters of the likelihood 
function were the number Nx of X particle sig- 
natures and their velocity vx- 

The result of such an analysis of the KAR- 
MENl data (1990-95), the negative logarithm of 
the likelihood —InL is shown in Figure ||a) as a 
function of one of the two free parameters, vx- 
A clear minimum at vx = 4.89m//is is evident. 
The characteristic side maxima result from the 
double pulse structure of pion decays as under- 
lined by Monte Carlo studies with high statistics 
(MC in Fig. |^) of X particle signals. Hence, the 
working hypothesis gave a consistent description 
of the anomaly in the time spectrum with a signal 
strength of Nx = 57 ± 25 excess events. 

One of the major objectives of KARMEN2 has 
been the investigation of this hypothesis, taking 
full advantage of the significantly reduced level of 
cosmic background. An initial KARMEN2 data 
sample sufficiently large for a detailed maximum 
likelihood analysis was collected from February 
1997 through February 1999. The corresponding 
maximum likelihood result —InL of these first 2 
years is given in Figure ^) having its minimum at 
Vx = 4.96 m/^s with a signal strength of Nx — 
17 ± 11. Scaling from the KARMENl analysis, 
one would have expected fq'^P'^'^^'^'^ — 22 ± 10. 
Though the minimum of the function —InL was 
not deep enough to form a statistically signifi- 
cant signal by KARMEN2 data alone, the fitted 
particle velocity and signal strength were in good 
agreement with the earlier KARMENl data, indi- 



cating a continuous accumulation of X signatures 
even under modified experimental systematics. 

However, an updated analysis of all KAR- 
MEN2 data collected so far (February 1997 
through March 2000) reveals no distortion of the 
time spectrum anymore. Despite an increase of 
pion and consequently neutrino production at the 
ISIS target by more than 50%, the number of ex- 
cess events in the time spectrum between 3.1 /zs 
and 4.1/xs after beam-on-target determined by 
subtraction of a fitted neutrino signal (exponen- 
tial decrease with r = 2.197^s) and constant 
background is now Nx = 11 ± 12.4. This allows 
to set an upper limit of Nx < 30 at a 90% con- 
fidence level. This has to be compared with an 
expected excess of ]\['^vected _ _|_ ^4 extrapo- 
lated from the KARMENl data. In addition to 
the time spectrum, the twodimensional maximum 
likelihood analysis shows no significant minimum 
of the negative logarithmic likelihood function at 
a velocity near vx — 4.96m//is (see Fig. ||c). Fix- 
ing the velocity vx at the above value results in a 
signal strength of Nx — 14 ± 12 events in contra- 
diction to the value of jsf'^P'^'^*'^'^ = 34 ± 15 extra- 
polated from the KARMENl likelihood analysis. 

To summarize, the time spectrum of the en- 
tire KARMEN2 data after 3 years of data tak- 
ing is in good agreement with the expectation 
of neutrino induced events with an exponential 
shape superimposed on a flat cosmic induced 
background. There is no distortion in the time 
spectrum nor any unexplained signal in a detailed 
maximum likelihood analysis as observed for the 
KARMENl data. As there was no change since 
1997 in data taking nor in the reduction and eval- 
uation of the data, the new results of the ana- 
lysis are due to the increased data amount and 
statistics only. On the other hand, the clear dis- 
tortion of the KARMENl spectrum remains un- 
changed and unexplained. While analysing the 
total KARMENl-h2 data set stiU yields a signi- 
ficant minimum in —InL at vx — 4.90m//is, the 
working hypothesis of n^—^ X has to be re- 

jected on the basis of the incompatibility of the 
'signal' strength seen by KARMENl and KAR- 
MEN2. 

KARMEN will continue to take data for an- 
other year, mainly to increase the sensitivity of 



the oscillation search, but also to monitor the 
time spectrum of single prong events. 
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Figure 5. Likelihood as a function of the velocity 
vx for KARMENl(a), initial KARMEN2(b) and 
all KARMEN2(c) data (see text for details). 



